Population redistribution in optically trapped polar molecules 
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Abstract. We investigate the rovibrational population redistribution of polar molecules in the electronic 
ground state induced by spontaneous emission and blackbody radiation. As a model system we use opti- 
cally trapped LiCs molecules formed by photoassociation in an ultracold two-species gas. The population 
dynamics of vibrational and rotational states is modeled using an ah-mitio electric dipole moment function 
and experimental potential energy curves. Comparison with the evolution of the ii"=3 electronic ground 
state yields good qualitative agreement. The analysis provides important input to assess applications of 
ultracold LiCs molecules in quantum simulation and ultracold chemistry. 
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Ultracold gases of dipolar molecules have been pro- 
posed as candidates for the exploration of quantum phases 
in dipolar gases [T] , the development of quantum compu- 
tation techniques [5], or precision measurements of fun- 
damental constants [3] . Recently tremendous progress has 
been made in the realization of cold and ultracold gases 
made out of polar molecules |4I5I6I7) . The presence of an 
electric dipole moment (EDM) in the electronic ground 
state leads to a coupling of the internal molecular state to 
the thermal environment via black-body radiation (BBR). 
Thus if the system is prepared in a certain internal state, 
this state will be altered after a finite time. In particular 
this yields a limited lifetime even for molecules prepared in 
the absolute rovibrational ground state. The thermaliza- 
tion of cold OH and OD molecules with the environment 
was observed experimentally [8] and theoretical studies 
were performed for a number of polar species J9il0j . For 
the molecular ions HD+ or MgH+ rotational cooling re- 
lying upon fast BBR-driven rotational relaxation was re- 
cently demonstrated by [TI] and [H], respectively. As a 
second consequence of the existence of an EDM, spon- 
taneous decay of electronic ground state molecules into 
more deeply bound rovibrational levels is possible. For the 
molecular ion HD"^ radiative decay leads to vibrational 
relaxation into w" = in less than 1 s [13114) . Radiative 
transitions rates for polar alkali dimers have been calcu- 
lated |15ll6j , but were generally considered to be too small 
to be of experimental relevance. 

Here we model the population redistribution in LiCs 
in the electronic ground state induced by spontaneous 
emission and BBR. We compare this simulation to mea- 
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surements on an ultracold gas of LiCs molecules in the 
electronic ground state, confined in an optical dipole trap 
insensitive to the internal state of the molecules. An ul- 
tracold gas of LiCs, the most polar alkali dimer [17j . in 
deeply bound levels of the X^ E'^ state can be efficiently 
formed by photoassociation |7I18) . We have measured the 
EDM of deeply bound ground state levels to be as high as 
5.5Debye I19|, making this molecule a promising candi- 
date for the realization of an ultracold gas with strong 
dipolar interactions. Here we observe significant dipole 
transition rates between vibrational ground state levels, 
enhanced by the strong dependence of the dipole moment 
function on the internuclear distance which is associated 
with its large average value. 



1 Ab-initio modeling of transition rates 

A transition between levels i and j with wave functions 
^i and ^j is proportional to the dipole matrix element 
pL^j = {^i\^{R)\\I/j)^ where pl{R) is the EDM function. For 
vibrational transitions only the variation of // with the 
internuclear distance it! leads to a finite transition dipole 
moment. The EDM function was calculated ab-initio and 
published previously [T7]- Rovibrational wave functions 
and overlap integrals are calculated with Level 8.0 [20] 
using a precise experimental potential energy curve [21) . 
The resulting EDMs of low lying vibrational levels are in 
very good agreement with measured values [191 . For mod- 
eling the spontaneous decay we calculate the Einstein A 
coefficients 



A,, 



A. 



Btt^S*, 



v' J' ,v" J" 



J', J" 



3£o^(2J' + l 



~3 2 

'^v'J',v"J"l^v'J',v"J" 



J. Deiglmayr^' et al.: Population redistribution in optically trapped polar molecules 



with the transition dipole moment fj.v'j',v"j" (in SI units), 
a transition energy i'v'j',v"j" (in wave numbers), and the 
Honl-London (HNL) Factor Sj'j» as defined in Ref. f^ . 
The coupling to BBR from the surrounding walls of the 
vacuum chamber is given by the Einstein B coefficients 
for stimulated absorption and emission multiplied by the 
spectral density of the BBR from an environment with 
given temperature T as 
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In Tab. [U the resulting rates at T == 293 K for the 
lowest vibrational levels are given. The given value for 
-B^'=i,i)"=o=^u'=o,i>"=i corresponds to a depopulation time 
constant of 65.4 s for molecules in the absolute ground 
state, in good agreement with previously calculated val- 
ues [9 10 . It is noteworthy that at cryogenic temperatures 
(77 K) this depopulation time constant is lengthened to 
1.3x10^ s. Already at room temperature rotational heat- 
ing due to BBR is greatly suppressed by the low spec- 
tral density at the relevant transition frequencies. Only 
for J" >12 rotational heating rates become comparable 
to experimental loss rates from the trap, which are on the 
order of tens of seconds. Rates for rotational transitions 
in the vibrational ground state are listed in Tab. [51 

The evolution of an initial distribution of populated 
ground state levels is finally given by the following rate- 
equation for the population Ni in a level i (the single index 
i enumerates rovibrational levels w", J" in energetic order- 
ing) 

+ E ^,.iVj(^) - E ^y^^O - ^^^W (1) 

where _r is a decay constant accounting for losses from the 
trap. The value for r=21s is taken from measurements 
(see below) and is assumed to be independent of v" . 

In the formation of ultracold molecules by photoasso- 
ciation (PA) [23124) , the initial population of ground state 
levels is determined by the Einstein A coefficients for the 
spontaneous decay from the excited PA state. For the B^TT 
state of LiCs we have calculated this initial population for 
different PA levels and found agreement with experimen- 
tal observations p^. The Frank-Condon factors between a 
given level in the B^TT state and levels in the X-'^Z'+ state 
show two maxima corresponding to a maximized overlap 
at the inner respectively the outer classical turning point 
of the vibrational motion. This leads to a typical bimodal 
distribution of the populated ground state levels with sig- 
nificant populations in deeply bound levels below v" < 8 
and in a small band of higher lying vibrational levels which 
shifts with the employed photoassociation resonance. The 
ground state population after PA via B^7T,ti'=26, J'=l is 
given in Tab. [T] for the deepest bound vibrational levels. 
The initial rotational distribution is determined by the 
HNL factors for the decay from the negative parity com- 
ponent of the excited PA level (due to the low collision 
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Fig. 1. Solution of Eq. [TJfor the population in the low- 
est nine vibrational levels of the X^ S^ after spontaneous 
decay from B^7T,w'=26,J'=l. The population is given as 
percentage of the initial total ground state population. 
Shown is also the population (rescaled to unity) after 
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e5 X A^l ~ 200 s. 



temperature of the atoms in the experiment s-wave scat- 
tering dominates iTl ) , as verified experimentally by deple- 
tion spectroscopy. Starting from this initial distribution 
the rate equation Eq. ([1]) is integrated numerically in or- 
der to derive the time evolution of the population, where 
convergence of the results is reached for time steps smaller 
than 50 ms. As can be seen in Fig. [1] the initial popu- 
lation distribution in the lowest vibrational levels of the 
X^Z'+ state is strongly modified. Most levels show either a 
steady increase or decrease essentially due to vibrational 
redistribution with adjacent states. As an exception the 
population in the v"=2i level exhibits a clear maximum 
after roughly 3 s, presenting a concise feature for testing 
the validity of the model. The time evolution of the pop- 
ulation in this level is shown in more detail in Fig. [5] a) 
together with the individual contributions from the rota- 
tional sublevels. The significant increase of the population 
in w"=3,J"=l, being caused by the decay of molecules 
from v" —A and the cycling of population between neigh- 
boring levels due to BBR with comparable rates, indicates 
stronger dynamics of the process than visible from the vi- 
brational populations only. 

We test the sensitivity of the model to the employed 
ab-initio dipole moment function by repeating the calcula- 
tion using a different dipole moment function published in 
Ref. [23 • The outcome of the simulation is modified only 
slightly and shows the same qualitative time evolution for 
the total population in vibrational levels v" < 8. 

As an alternative test we replace the first order loss 
term rNi{t) in Eq.[T]by a second order loss term where the 
values for the rate constant /3LiCs-LiCs=2xlO^^" cm^/s, 
the initial mean molecular density no,LiCs=2xlO* cm"^. 
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2 


3 
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5 


6 


7 


8 


Pop. %] 


0.0 


0.0 


0.5 


2.4 


5.1 


2.9 


0.0 


2.5 


0.3 


A„„+i,„„ 10-Vs 


22.7 


41.8 


57.4 


69.6 


78.6 


84.5 


87.5 


87.9 


85.8 


B„"+i,„" lO^Vs 


15.7 


29.3 


41. 


50.6 


58.1 


63.6 


67.1 


68.6 


68.3 


B„„ [cm-'] 


0.0 


182.7 


363.4 


542.0 


718.6 


893.2 


1065.7 


1236.1 


1404.4 



Table 1. Relative population strength of the relevant lowest vibrational levels of the X^Z'+ state after spontaneous 
decay from B^77,w'=26; total vibrational rates ^i,'.,," and Byi,yii (for T = 293 K), summed over rotational quantum 
numbers in A^iji^yiijn and Byij/jjujn] energy of vibrational levels Eyii. Rates for transitions with Aw" > 2 are 
significantly smaller and are thus not explicitly given, but are included in the simulation. 
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1.6 


5.9 


14.6 


29.2 


51.2 


82.2 


123.8 


177.3 
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0.7 


1.5 


2.5 


3.9 


5.5 


7.3 


9.5 


11.9 



Table 2. Rate constants Ay'j'yj" and i3u',,/',i,",,/" (for T — 293 K) 
ground state. The rotational constant for v"=0 is 0.1874 cm~^. 



for pure rotational transitions in the vibrational 



and the initial total number of molecules A^o.tot^SxlO'^ 
are taken from experimental estimates (see below): 

+ ^B,,7V,(i)-^B,,iV,(t) (2) 
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As shown in Fig. [5] b) this leads to a significant faster 
decay which masks the redistribution process in the inte- 
grated vibrational signal while the dynamics remain visi- 
ble in the single rotational components. 



2 Experimental observation of population 
dynamics 

In the experiment ultracold LiCs molecules are formed 
by photoassociation (PA) |23I24) of ^^^Cs and "^Li atoms 
trapped simultaneously in a quasi electrostatic trap (QUEST) 
After the removal of remaining atoms a pure molecular 
sample can be kept in the trap for variable times. The 
QUEST is formed by the single focus of a CO2 laser (power 
64(4) W, waist wo=89(4) [im) emitting at 10.6 \xvn. De- 
tails of the QUEST setup can be found in [26] and refer- 
ences therein. The resulting trap depths and frequencies 
for the different species are given in Tab. [31 They are cal- 
culated from measured values for beam power and waist, 
experimental values for atomic polarizabilities [37], and 
ab-initio calculated values for the molecular polarizabili- 
ties [5S]. 

Despite its significant detuning, the optical trapping 
field could drive Raman transitions between neighboring 
levels due to its high intensity. The rate for such transi- 
tions can be estimated from a simple classical model of the 
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Fig. 2. Simulated time evolution of the population in 
X^Z'+,w"=3. Shown are the total signal (solid line) and 
the population of the rotational sublevels (dashed lines, 
labeled by J"), a) Simulation based on a first order loss 
mechanism (eq. [Ij , b) simulation based on a second order 
loss mechanism (eq. [5]). 
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Species 


Li 


Cs 


LiCs(x'i;+,i)"=o) 


LiCs(x^i:+y=42) 


LiCs(a^X'+,i;"=0) 


Trap depth [/iK] 

^rad/2n [l/s] 
Ula:./2n [l/s] 


194(21) 
1729(166) 

47(7) 


479(53) 
623(60) 

17(2) 


435(48) 

579(55) 

16(2) 


700(77) 

735(70) 

20(3) 


744(82) 

758(73) 

20(3) 



Table 3. Trap depths and frequencies of the QUEST for the different species. The given error is based on the exper- 
imental uncertainties only, no error on the theoretical values was assumed. X^Z'+,ti"=42 experiences the maximum 
trap depth of all X^S^ levels. 



Raman process [33]. Using the ab-initio calculated polar- 
izability tensor for the X^E'^ state of LiCs [IS] and the ex- 
perimental laser intensity we derive transition rates which 
are several orders of magnitude below the transition rates 
for spontaneous decay and BBR transitions. Thus vibra- 
tional or rotational heating by the trapping CO2 laser can 
be neglected even on longer timescales. 

For PA of LiCs molecules in the dipole trap we load 
first Cs atoms from a Zeeman slowed atomic beam into 
a conventional MOT setup. After a red-detuned molasses 
phase 8x10^ Cs atoms are transferred into the QUEST at 
a density of 1 x 10^^ cm^'^ and a temperature of 33(4) /iK. 
Afterwards the Li MOT is loaded also from a Zeeman 
slowed atomic beam. After magnetic compression typically 
2x10^ Li atoms are transferred into the dipole trap at a 
density of 5 x 10^ cm~^ . The comparatively low density has 
to be attributed to the relatively small initial number of 
2x10^ atoms in the MOT, and to a low transfer efficiency 
due to the small phase-space overlap between MOT and 
QUEST. Li and Cs atoms are then kept together in the 
trap for 500 ms to reduce the Li temperature by thermal- 
ization with the Cs ensemble {^l . The background gas in 
the science chamber at a pressure of around 2 x 10~^^ mbar 
leads to atomic lifetimes on the order of 30 seconds. 
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Li(2i 5,n) 
+Cs(6pXj 



+Cs(6/5,„) 



Fig. 3. Exemplary processes in the formation and detec- 
tion of a ground state LiCs molecule: a) photoassociation 
from the continuum into the B^il state, b) spontaneous 
decay into the X.^E~^ , c) two-photon ionization via inter- 
mediate levels in the B^TJ state. Redistribution of popula- 
tion: d) spontaneous decay, e) transitions driven by black- 
body radiation (for clarity only transitions with Av" = ±1 
are shown). 



After these preparatory steps molecules are formed 
and detected as illustrated in Fig. [31 First PA is per- 
formed by a light pulse (300 ms long) from a tunable cw 
TiSa laser (typically 500 mW, focused to 100 fim). The 
excited LiCs molecules decay spontaneously within few 
tens of nanoseconds into different electronic ground state 
levels. Short pulses of resonant light are used to remove 
remaining atoms, then all laser beams are blocked by me- 
chanical shutters. After a variable hold time in the dipole 
trap, molecules in given vibrational ground state levels 
are ionized by resonant-enhanced multi-photon ionization 
(REMPI) and the resulting ions are detected in a time- 
of-flight mass spectrometer. From the number of detected 
ions we deduce an overall number of roughly 5x10'^ molecules21(l) seconds (comparable to the lifetime of pure atomic 



peak density of nucs ~6x 10^ cm~'^. For molecules in high 
vibrational levels the density could be up to 60% higher 
due to the increased polarizability. 

The slowest decay from the trap is observed for molecules 
in high lying levels, shown in Fig. 21 Assuming a first order 
decay process a fit to the data yields a time constant of 



in a broad distribution of internal states [18] which are ini- 
tially stored in the trap. The translational temperature of 
the molecules is mostly determined by the temperature 
of the heavier Cesium and can thus be estimated to be 
around 50/xK. The density of the molecular ensemble is not 
easily defined, as the trapping potential depends on the vi- 
brational state. However deeply bound LiCs molecules ex- 
perience a similar trapping potential as Cs atoms and one 
can thus derive the molecular density from the measured 
atomic density as nucs = A^LiCs/A^Cs "-Cs- This yields a 



samples in the trap). At our estimated molecular density 
bi- molecular collisions might also become significant. Fit- 
ting the decay curve with a pure two-body decay yields a 
rate coefficient of /3LiCs-LiCs=2xl0~^" cm^/s for inelastic 
LiCs-LiCs collisions (making the above discussed assump- 
tions about the molecular density). This is a realistic value 
when compared with rate coefficients found in other exper- 
iments. For collisions of KRb molecules in deeply bound, 
distinguishable states a rate coefficient of 2 x f O^^'' cm'^/s 
[31j was measured. For RbCs molecules in different weakly 
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Fig. 4. Ion yield at 16893. 3cm~^ and high intensity (■~350 
J/cm^) after a varying hold time in the QUEST. Under 
these REMPI conditions molecules in levels X^U+,v">30 
can be detected. Every point is an average of nine mea- 
surements, shown uncertainties are the standard deviation 
of the measurement. The solid (dashed) line shows a fit of 
a first-order (second-order) decay process to the data. 
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Fig. 5. Ion yield from v"=3 molecules ionized at 
16886. 3cm~-'^ and low intensity (~2 J/cm^)after a vary- 
ing hold time in the QUEST. Every point is an average of 
15 measurements, given uncertainties are the standard de- 
viation of the mean value. The dotted (dashed) line shows 
the simulated signal for the adjusted first (second) order 
model as discussed in the text. 



bound triplet states a value on the order of 1 x 10^^" cm^/s 
was found [22] ■ However our data doesn't allow us to dis- 
tinguish between a first or second order decay process. As 
the observed loss rate is comparable to the atomic ones, 
we cannot exclude that the observed molecular trap loss 
is due to collisions with background gas or another yet 
unidentified loss mechanism. 

At typical REMPI intensities of 2 J/cm^ the observed 
REMPI resonances have a width of 2-3 cm^^ (full width at 
half maximum), which is very likely due to power broad- 
ening of the REMPI transitions (the line width of the 
ionization laser is 0.07 cm~^, the spread of a typical ini- 
tial rotational distribution ~0.5 cm~^). Simulations of the 
REMPI line profile after different storage times (and thus 
different rotational distributions) show that the ionization 
signal at the initial peak position of the REMPI signal is 
a reliable measurement for the total population in a vi- 
brational level. 

Fig. [5] shows the ion yield from molecules in the vibra- 
tional level X^Z'+,d"=3. We used depletion spectroscopy 
to check that at the employed REMPI energy of 16886.3 cm 
the detected ions originate indeed predominantly from this 
ground state level. The experimental curve shows the same 
qualitative behavior as seen for the simulation assuming a 
first order loss process (Eq.[T]): first an increase is observed, 
then a maximum is reached after around 3s followed by a 
decay of the signal. However the modulation of the exper- 
imental curve is significantly stronger than in the simula- 
tion. The fast decay for hold times longer than 3s is in fact 
better described by the model assuming a second order 
loss process. In order to adjust the model to the data we 
vary only the relative initial population in the vibrational 
levels v"=3 and v"—4: which determines most strongly the 
modulation of the population in v"—S. The initial popu- 
lation is the most sensitive quantity in the model as it 



is mostly determined by the Franck-Condon (EC) factors 
between the PA level and the electronic ground state lev- 
els which are sensitive to small deviations of the potential 
energy curves. The subsequent dynamic of the population 
in the ground state on the other hand is dominated by 
transitions between neighboring levels in the X^Z""*" state 
for which the rates are less sensitive to the exact shape 
of the potential. Fig. [5] shows the result of simulations 
where the ratio Ny"=4/Ny"=3 was adjusted to reproduce 
the experimental modulation strength. For the first order 
model of Eq. [T]the ratio was increased by a factor of 2.7, 
for the second order model of Eq. [2] by a factor of 5.0. 
The latter simulation reproduces the initial modulation of 
the experimental data very well. However the necessary 
adjustment of the starting conditions is stronger than one 
would expect for EC factor calculation from accurate ex- 
perimental potentials. The experimentally observed decay 
of the population in u"=3 is still faster than expected 
from our measurements of the decay for high lying vibra- 
tional levels. This could be an indication for an increase 
of the cross section for inelastic LiCs-LiCs collisions with 
-^increasing binding energies as predicted for the case of 
RbCs [33,. 



3 Conclusion and outlook 

We have observed population redistribution caused by BBR 
and spontaneous decay for ultracold LiCs molecules stored 
in an optical dipole trap. Similar rate calculations for other 
polar molecules which have been formed in the ultracold 
regime, KRb 6^ and RbCs 34,, show that the timescales 
for redistribution processes in these systems is orders of 
magnitude larger ^9ll0) . Due to the large electric dipole 
moment of ground state LiCs molecules and the associ- 
ated strong variation of the dipole moment function with 
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the internuclear distance, radiative and BBR-driven dipole 
transitions occur in this system on a timescale relevant for 
this and future experiments. In the specific case of LiCs, 
molecules prepared in the rovibrational ground state will 
be pumped out of this level by BBR-driven transitions 
with a time constant of 65 s. If in future experiments 
this imposes a limitation, this value could be significantly 
lengthened by cooling the apparatus. The comparison of 
our experimental results with models based on ab initio 
calculated rates may indicate the presence of a bimolecu- 
lar loss with a rate coefficient exceeding 2x10"^" cm^/s 
for deeply bound LiCs molecules. 

This work is supported by the DFG under WE2661/6-1 in the 
framework of the Collaborative Research Project QuDipMol 
within the ESP EUROCORES EuroQUAM program. JD ac- 
knowledges partial support of the French- German University. 
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